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A  numerical  study  using  crystal  plasticity  finite  element  method  was  performed  in  order  to  investigate  the 
influence  of  the  grain  boundary  morphology  (stair-case  or  smooth/flat)  and  the  selection  of  different  mesh 
types  (hexahedral  and  tetrahedral  meshes)  on  the  distribution  of  simulated  plasticity  heterogeneities.  A 
bicrystal  (a  hard  grain  and  a  soft  grain)  sphere  embedded  in  a  cubic  grain  was  utilized  to  perform  the 
current  numerical  study.  The  volume  averaged  responses  for  the  simulated  plasticity  values  (von  Mises 
stress,  total  accumulated  slip  system  shear  and  maximum  accumulated  slip  system  shear)  showed  no 
significant  dependence  on  the  grain  boundary  representation  and  mesh  type.  However,  the  bicrystal  rep¬ 
resented  by  uniform  hexahedral  meshes  with  stair-case  grain  boundaries  tends  to  show  more  extreme 
local  plasticity  values  (in  the  head  and  tail  of  the  distribution),  particularly  populated  in  the  grain  bound¬ 
ary  region.  Differences  in  local  plasticity  values  between  the  stair-case  and  the  smooth  (or  flat)  grain 
boundary  cases  and  between  the  hexahedral  and  the  tetrahedral  mesh  cases  were  largest  in  the  grain 
boundary  region,  and  the  degree  of  the  differences  depended  upon  the  plasticity  susceptibility  of  the 
grain.  The  simulation  results  suggest  that  the  grain  boundary  represented  by  the  stair-case  morphology 
can  be  a  source  and/or  a  sink  for  local  extreme  plasticity,  compared  to  the  grain  boundary  represented 
by  the  smooth  (or  flat)  morphology. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Crystal  plasticity  finite  element  method  (CP-FEM)  approaches 
have  been  widely  used  to  model  elasto-viscoplastic  or  plas¬ 
tic  responses  of  polyctystalline  metals  at  the  mesoscopic  scale, 
wherein  deformation  heterogeneities  due  to  anisotropic  responses 
of  differently  oriented  grains  are  of  main  interest.  In  general, 
such  mesoscopic  CP-FEM  approaches  utilize  3D  polycrystalline 
microstructures,  which  are  synthetically  generated  using  numer¬ 
ical  tools  [1-6],  experimentally  obtained  using  destructive  or 
non-destructive  techniques  [7,8],  or  semi-empirically  recon¬ 
structed  from  experimental  2D  images  [9-14].  Simulation  results 
are  usually  analyzed  in  order  to  clarify  the  spatial  distribution 
of  stress  (or  strain)  heterogeneities  (or  hot  spots)  [9,14,15],  the 
statistics  of  extreme  values  of  simulated  plasticity  [16]  and  its 
relation  to  microstructural  features  such  as  grain  boundaries  and 
precipitates  (if  present).  Here,  the  heterogeneity  of  plasticity  near 
grain  boundaries  draws  particular  attention  since  such  regions 
can  be  susceptible  to  plasticity,  depending  upon  the  level  of  the 
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deformation  compatibility  between  the  grains  neighboring  each 
other  [14,17]. 

In  reality,  grain  boundaries  in  polycrystalline  metals  are 
smooth  with  curvatures,  based  upon  the  energy  balance  from 
crystallography  and  chemistry  between  neighboring  grains.  3D 
microstructures  used  for  CP-FEM  simulations,  however,  are  mostly 
represented  in  uniform  cubic  meshes  (or  “voxels”),  which  causes 
grain  morphologies  inherently  represented  by  “stair-case"  bound¬ 
aries.  It  is  generally  expected  that  simulated  local  plasticity 
responses  can  be  different  between  3D  microstructures  having 
stair-case  grain  boundaries  and  smooth  (or  flat)  grain  boundaries, 
particularly  in  the  vicinity  of  grain  boundaries.  Furthermore,  even 
for  a  3D  microstructure  properly  meshed  to  have  smooth  (or  flat) 
grain  boundaries,  the  choice  of  a  mesh  type,  such  as  hexahedral  and 
tetrahedral  meshes,  may  also  influence  the  distribution  of  simu¬ 
lated  local  plasticity  responses.  It  is  surprising  to  know  that  only 
a  few  numerical  studies  have  been  actually  conducted  to  com¬ 
pare  crystal  plasticity  responses  between  different  grain  boundary 
representations  [4,5,18]  and  mesh  types  [19,20].  Even  for  those 
numerical  studies,  the  comparison  was  limited  to  the  global  aver¬ 
age  response  [5,18]  and  just  qualitative  local  responses  [4]. 

Increasing  the  mesh  density  can  visually  minimize  stair-case 
boundaries  of  grains  and  improve  the  discrepancy  in  local  plasticity 
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Fig.  1.  (a)  Schematic  illustration  of  the  tricrystal  geometry  used  in  the  present  study,  and  3D  cross-sectional  views  of  the  tricrystal  differently  represented  by  (b)  hex-meshed 
stepped,  (c)  hex-meshed  smooth  and  (d)  tet-meshed  smooth  geometries,  respectively. 


responses  between  the  smooth  and  the  stepped  cases  [21]  at 
a  certain  resolution  window.  However,  this  approach  can  have 
a  prohibitive  computational  expense  (in  computation  time  and 
resources),  and  still  does  not  address  the  fundamental  issue  associ¬ 
ated  with  the  difference  in  deformation  compatibility  and  resulting 
heterogeneity  between  the  smooth  (or  flat)  and  the  stair-case 
boundaries.  This  issue  and  associated  risks  need  to  be  clearly 
assessed  and  understood  when  the  CP-FEM  simulation  is  intended 
to  study  the  statistics  of  hot  spots  and  extreme  values  of  plastic¬ 
ity.  In  the  present  study,  preliminary  CP-FEM  simulations  were 
performed  using  a  simple  tricrystal  geometry  (a  bicrystal  sphere 
embedded  in  a  cubic  grain)  in  order  to  gain  basic  understand¬ 
ing  about  differences  in  simulated  mesoscopic  plasticity  responses 
and  heterogeneities  between  3D  microstructures  having  differently 
represented  grain  boundary  morphologies  and  mesh  types. 


2.  Simulation  details 

A  tricrystal  geometry,  which  is  a  bicystal  sphere  embedded  in 
a  cubic  grain,  was  generated  for  the  current  parametric  CP-FEM 
study.  Fig.  1(a)  schematically  illustrates  the  tricrystal  geometry. 
A  diagonally  flat  grain  boundary  was  introduced  to  the  bicrys¬ 
tal  inside  the  cubic  grain.  The  ratio  between  the  diameter  of  the 
bicrystal  sphere  and  the  edge  length  of  the  cubic  grain  was  main¬ 
tained  to  be  2:3.  The  spherical  bicrystal  separated  by  the  flat 


diagonal  boundary  is  intended  to  investigate  the  difference  in  local 
plasticity  responses  when  both  curved  and  flat  grain  boundaries 
are  represented  by  different  topologies  (smooth  and  stair-case 
topologies)  and  mesh  types  (hexahedral  and  tetrahedral  elements). 
Fig.  l(b)-(d)  shows  3D  cross-sectional  views  of  the  tricrystal  dif¬ 
ferently  represented  by  (b)  hex-meshed  stepped,  (d)  hex-meshed 
smooth  and  (d)  tet-meshed  smooth  geometries,  respectively.  For 
the  generation  of  the  stair-case  geometry  of  the  tricrystal  shown 
in  Fig.  1(b),  the  original  geometry  of  Fig.  1(a)  was  superimposed 
on  a  36  X  36  X  36  hexahedral  (uniform  cubic)  grid,  and  each  hex- 
element  was  assigned  to  a  corresponding  grain,  based  upon  its 
proximity  to  the  grain.  For  the  generation  of  smooth  geometries 
of  the  tricrystal  shown  in  Fig.  1(c)  and  (d),  all  three  grains  were 
meshed  independently  (with  hexahedral  and  tetrahedral  elements, 
respectively)  while  maintaining  the  smooth  and  diagonally  flat 
boundaries,  and  the  continuity  of  meshes  across  the  boundary. 

Fig.  2  shows  the  distribution  of  minimum  dihedral  angles  (^min) 
for  3D  mesh  elements  that  constitute  the  bicrystal  sphere  repre¬ 
sented  by  hex-meshed  smooth  (Fig.  1(c))  and  tet-meshed  smooth 
(Fig.  1(d))  geometries,  respectively.  Note  that  the  dihedral  angle  for 
the  perfect  equilateral  tetrahedron  is  70.5°.  In  Fig.  2,  the  average 
0niin  for  tet-  and  hex-mesh  elements  is  about  50°  and  80°,  respec¬ 
tively,  and  the  smallest  ate  about  20°  for  tet-mesh  elements 
and  about  58°  for  hex-mesh  elements.  It  indicates  that  quality  hex- 
and  tet-meshes  were  generated  for  smooth  bicrystal  geometries 
such  that  one  can  expect  minimized  numerical  errors  associated 
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Fig.  2.  Distribution  of  minimum  dihedral  angles  for  the  bicrystal  sphere  represented 
by  hex-meshed  smooth  (Fig.  1(c))  and  tet-meshed  smooth  (Fig.  1(d))  geometries, 
respectively. 


with  poor  quality  mesh  elements  [22].  Also,  an  extra  effort  was 
made  to  mesh  all  cases  in  Fig.  l(b)-(d)  to  have  a  computationally 
comparable  number  of  mesh  elements  such  that  the  computational 
cost  is  of  the  same  order  for  all  cases.  Total  numbers  of  elements 
and  corresponding  integration  points  for  each  case  are  summarized 
in  Table  1  for  each  grain.  Here,  linear  elements  (C3D8  for  hexahe- 
dral  and  C3D4  for  tetrahedral  elements  [22])  were  used  for  CP-FEM 
simulations  for  all  cases. 

The  elasto-viscoplastic  constitutive  model  for  FCC  crystals  was 
used  in  the  present  CP-FEM  simulations.  In  particular,  anisotropic 
elasto-viscoplasticity  was  applied  to  the  bicrystal  sphere  (grains 
2  and  3  in  Fig.  1(a)),  while  isotropic  elasto-viscoplasticity  was 
used  for  the  surrounding  cubic  grain  (grain  1)  in  Fig.  1(a).  For  the 
anisotropic  elasto-viscoplasticity  constitutive  model  for  grains  2 
and  3  in  Fig.  1(a),  the  power  law  and  Voce-type  parabolic  harden¬ 
ing  formulations  were  used  to  describe  the  slip  system  shear  rate 
()/“)  and  the  evolution  of  slip  resistance  (g): 


l/m 


sign(T^) 


and 


l  =  eol^V| 

gs  -  go 


(1) 

(2) 


where,  yo ,  and  m  are  the  reference  shear  rate,  the  resolved  shear 
stress  on  the  slip  system  a  and  the  strain  rate  sensitivity  parame¬ 
ter,  respectively.  In  Eq.  (2),  0o,  go  and  gs  are  the  initial  hardening 
rate,  initial  and  final  values  of  the  slip  resistance,  respectively.  Val¬ 
ues  of  these  parameters  used  in  the  present  simulations  are  listed 
in  Table  2.  Here,  Cn,  Q2  and  C44  are  elastic  constants.  Elastic  con¬ 
stants  in  Table  2  were  taken  from  those  for  pure  nickel.  However, 
all  other  input  parameter  values  in  Table  2  was  arbitrarily  chosen 


Table  1 

Total  numbers  of  elements  and  integration  points  (in  brackets)  for  the  triciystal. 


Hex- 

meshed 

stepped 

Hex- 

meshed 

smooth 

Tet- 

meshed 

smooth 

Grain  2 

3.6  (28.8) 

3.7  (29.8) 

29.5  (29.5) 

Grain  3 

3.6  (28.8) 

3.7  (29.8) 

29.4  (29.4) 

Grain  1 

39.4(315.6) 

38.4(307.3) 

321.1  (321.1) 

In  thousands. 


Table  2 

Values  of  input  parameters  used  for  anisotropic  elasto-viscoplasticity  (grains  2  and 
3). 


Vo 

0.001  s-> 

m 

0.03 

go 

200  MPa 

gs 

400  MPa 

Bo 

500  MPa 

Cn 

247  GPa 

Cl2 

147  GPa 

C44 

125  GPa 

to  represent  a  FCC  ciystal.  For  isotropic  elasto-viscoplasticity  for 
grain  1  in  Fig.  1(a),  values  for  Cn,  C12  and  C44  were  given  as  344, 
147  and  98GPa,  respectively,  such  that  the  elastic  anisotropy  fac¬ 
tor  (A  =  C44/0.5(Cii  -  C12))  becomes  1  for  isotropic  elasticity.  Upon 
plasticity,  slip  rate  (/“)  was  set  to  be  equal  for  all  twelve  slip  sys¬ 
tems  in  order  to  achieve  isotropic  plasticity  -  strictly  speaking, 
isotropic  plasticity  used  in  the  present  model  is  limited  to  isotropic 
plastic  flow  for  twelve  slip  systems  of  a  FCC  crystal. 

Constitutive  models  described  above  were  implemented  in  the 
commercial  finite  element  package  ABAQUS®  through  the  User 
MATerial  subroutine  UMAT  [22].  Simple  tension  (along  the  vertical 
direction  in  Fig.  1  (a))  boundary  conditions  were  applied  to  all  simu¬ 
lation  cases  given  in  Fig.  1(b)  and  (c),  with  the  strain  rate,  0.0001  s“' , 
up  to  the  target  axial  strain,  0.6%.  Here,  grain  2  and  grain  3  in  Fig.  1(a) 
are  oriented  such  that  their  [ill]  and  [123]  crystallographic  direc¬ 
tions  are  parallel  to  the  tensile  direction,  respectively.  This  sets 
grains  2  and  3  as  plastically  hard  and  soft,  respectively  (when 
pulled  along  the  tensile  direction).  Results  from  all  simulation  cases 
(Fig.  l(b)-(d)),  particularly  for  grains  2  and  3,  were  quantitatively 
compared  by  their  volume  averaged  global  responses,  the  statis¬ 
tical  distribution,  and  upper  and  lower  plasticity  values  (viz.,  hot 
and  cold  spots).  Here,  three  selected  parameters,  the  von  Mises 
stress  (ctm).  the  total  accumulated  plastic  shear  strain  (Ey“)  from 
all  slip  systems,  and  the  maximum  slip  shear  strain  accumulated  in 
a  slip  system  (Kmax)  were  chosen  to  compare  local  plasticity  values 
between  the  cases.  For  the  calculation  of  E/“  the  increment  of  slip 
system  shear  (Ay")  was  summed  over  all  twelve  slip  systems  at 
each  simulation  time  step  (|EAy“|4f),  and  Ey“  was  determined 
by  the  summation  of  all  |EAy“|At’s. 

3.  Simulation  results  and  discussion 

Fig.  3(a)-(c)  shows  responses  of  CTm,  Ey“  and  values 
averaged  over  grains  2  and  3,  respectively,  as  a  function  of  the 
global  axial  strain  for  hex-meshed  stepped,  smooth  and  tet-meshed 
smooth  cases.  In  Fig.  3(a)  CTm  values  averaged  over  all  three  grains 
were  also  plotted  as  a  function  of  the  global  axial  strain  for  com¬ 
parison.  The  simulation  result  shows  no  significant  differences  in 
grain-averaged  ctm.  and  y"^^  responses  as  a  function  of  the 
global  axial  strain  for  all  simulation  cases.  In  Fig.  3(a),  the  <T|vi 
response  averaged  over  all  three  grains  shows  no  sensitivity  to 
the  choice  of  stair-case  or  smooth  grain  boundaries  and  hexahe- 
dral  or  tetrahedral  meshes.  The  maximum  percentage  difference 
in  grain-averaged  crM.  2/“  and  y"^^  values  at  0.6%  axial  strain 
was  0.4%,  3%  and  3.3%  for  grain  2  and  1%,  0.7%  and  1%  for  grain  3, 
respectively.  The  results  clearly  show  that  the  average  responses  of 
simulated  o'M,  Ey“  and  values  are  not  significantly  influenced 
by  the  selection  of  the  grain  boundary  representation  (stair-case 
or  smooth  boundaries)  or  the  mesh  type  (hexahedral  or  tetra¬ 
hedral  meshes).  The  statistical  distribution  of  all  simulated  o'm, 
Ey“  and  values  (taken  from  all  integration  points)  at  0.6% 
strain  for  grains  2  and  3  is  plotted  in  a  semi-log  scale  in  Fig.  4  for 
hex-meshed  stepped,  hex-meshed  smooth  and  tet-meshed  smooth 
cases.  In  Fig.  4,  the  statistical  distribution  of  crM.  S/“,  values  is 
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Fig.  3.  Simulated  volume  averaged  responses  of  (a)  von  Mises  stress,  (b)  total  accumulated  shear  strain,  and  (c)  maximum  slip  system  shear  strain  as  a  function  of  global 
axial  strain. 
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Fig.  4.  Statistical  distribution  of  all  simulated  ctm.  >^3^  values  (taken  from  all  integration  points),  as  semi-log  plots,  for  grain  2  ((a)-(c))  and  for  grain  3  ((d)-(f))  for 
hex-meshed  stepped,  hex-meshed  smooth  and  tet-meshed  smooth  cases  at  0.6%  strain. 
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Fig.  5.  All  extreme  values  of  (a)  ctm.  (b)  S)/“  and  (c)  (from  Fig.  4)  plotted  as  a  function  of  the  normalized  distance  from  the  grain  boundary  for  the  hex-meshed  stepped 
grains  2  and  3.  The  normalized  distance  from  the  grain  boundary  was  determined  by  normalizing  the  calculated  grain  boundary  distances  with  the  maximum  grain  boundary 
distance. 


very  similar  between  the  hex-meshed  smooth  (green  squares)  and 
the  tet-meshed  smooth  (red  triangles)  cases.  However,  the  hex- 
meshed  stepped  case  (blue  squares)  tends  to  show  more  excessive 
O'M.  and  values,  both  in  the  tail  and  head  of  the  dis¬ 
tribution,  compared  to  the  hex-meshed  and  tet-meshed  smooth 
cases.  These  excessive  cold  and  hot  (extreme)  spots  in  the  hex- 
meshed  stepped  case  are  particularly  apparent  for  the  CTm. 
and  distribution  in  grain  2  (Fig.  4(a)-(c)),  and  the  CTm  distri¬ 
bution  in  grain  3  (Fig.  4(d)).  The  E)/“  and  distribution  in  grain 
3  (Fig.  4(e)-(f))  is  almost  identical  for  the  hex-meshed  stepped, 
hex-meshed  smooth  and  tet-meshed  smooth  cases. 

The  results  from  Fig.  4  indicate  that  a  grain  represented  by  stair¬ 
case  grain  boundaries  can  cause  excessively  low  and/or  high  local 
plasticity  values,  compared  to  a  grain  represented  by  smooth  (or 
fiat)  grain  boundaries.  The  presence  of  these  excessive  extreme 
values  depends  upon  the  plastic  susceptibility  of  the  grain  and 
the  choice  of  a  plasticity  parameter.  In  order  to  clarify  how  these 
extreme  plasticity  values  are  spatially  distributed  in  the  hex- 
meshed  stepped  grains,  all  extreme  CTm.  S)/“  and  values, 
which  are  smaller  (or  larger)  than  the  minimum  (or  maximum)  val¬ 
ues  for  the  hex-meshed  and  tet-meshed  smooth  cases,  were  plotted 
as  a  function  of  the  normalized  distance  from  the  grain  boundary 
in  Fig.  5.  Here,  volume  fractions  of  these  excessive  upper  and  lower 
extreme  values  found  in  the  hex-meshed  stepped  grains  were  less 
than  1%.  It  is  clear  in  Fig.  5  that  all  upper  extreme  ctm.  and 
values  in  the  hex-meshed  stepped  grains  are  distributed  very  near 
to  the  grain  boundary,  as  are  most  of  the  lower  extreme  values. 

It  seems  to  be  apparent  from  Figs.  4  and  5  that  the 
hex-meshed  stepped  grain  geometry  tends  to  have  upper 
and  lower  extreme  spots  of  selected  parameters  that  exceed 
those  from  the  hex-meshed  and  tet-meshed  smooth  grain 
geometries.  Also,  these  extreme  spots  were  mainly  found  near 
the  grain  boundaries.  This  implies  that  local  plasticity  values  can 
be  underestimated  or  overestimated  in  the  grain  boundary  region 
when  the  grain  is  represented  by  stair-case  boundaries,  in  compar¬ 
ison  to  the  cases  that  the  grain  is  represented  by  smooth  (or  flat) 
boundaries.  The  degree  of  the  underestimation  and  overestimation 
seems  to  be  dependent  upon  how  the  deformation  compatibility 
interacts  between  neighboring  grains  in  the  stair-case  boundary 
region.  Here,  a  quantitative  comparison  of  local  plasticity  values 
between  the  cases  is  expected  to  provide  the  useful  information. 
However,  the  exact  point-to-point  quantitative  comparison  in  plas¬ 
ticity  parameter  values  between  the  stepped  and  the  smooth  grain 
geometries  is  not  feasible  because  local  grain  boundary  regions 
occupied  by  the  stepped  grain  and  the  smooth  grain  are  different. 
To  roughly  estimate  the  local  plasticity  parameter  value  differences 
between  the  cases,  the  nearest  point-based  selected-parameter 
value  comparison  was  performed  for  the  hex-meshed  stepped  and 


hex-meshed  smooth  grain  geometries,  and  for  the  hex-meshed 
smooth  and  tet-meshed  smooth  grain  geometries.  In  Fig.  6  the 
average  percentage  differences  in  local  selected-parameter  val¬ 
ues  between  the  hex-meshed  stepped  and  hex-meshed  smooth 
cases  (Fig.  6(a)-(c)),  and  between  the  hex-meshed  smooth  and 
the  tet-meshed  smooth  cases  (Fig.  6d-f),  are  plotted  as  a  func¬ 
tion  of  normalized  distance  from  the  grain  boundary  for  grain  2 
(closed  symbols)  and  grain  3  (open  symbols)  at  0.6%  strain.  In 
Fig.  6(a)-(c),  the  figure  legend  “Hex,  Stepped  >  Hex,  Smooth”  (or 
“Hex,  Stepped  <  Hex,  Smooth”)  indicates  the  average  percentage 
difference  profiles  (plotted  as  a  function  of  normalized  distance 
from  the  grain  boundary)  for  all  local  o'M,  and  val¬ 
ues  whose  magnitudes  are  larger  (or  smaller)  in  the  hex-meshed 
stepped  geometry  than  in  the  hex-meshed  smooth  geometry.  The 
same  notation  applies  to  the  figure  legend  “Hex,  Smooth  >Tet, 
Smooth”  and  “Hex,  Smooth  <Tet,  Smooth”  in  Fig.  6(d)-(f).  Here, 
several  things  need  to  be  pointed  out  from  Fig.  6. 

First,  the  average  percentage  differences  in  local  plasticity  val¬ 
ues  (hereafter,  for  simplicity  these  will  be  referred  to  as  Aom. 
AEy“  and  Ay"^^  for  average  percentage  differences  in  local  ctm. 
Ey“  and  values  between  the  cases,  respectively)  tend  to  be 
larger  in  the  grain  boundary  region,  compared  to  those  in  the 
grain  interior  region,  for  all  comparison  cases.  This  trend  is  under¬ 
standable  when  comparing  local  plasticity  differences  between 
the  hex-meshed  stepped  and  the  hex-meshed  smooth  cases,  since 
the  mesh  morphological  difference  is  greatest  only  along  the 
grain  boundary.  Perhaps,  the  different  deformation  incompatibil¬ 
ity  associated  with  different  grain  boundary  morphologies  leads  to 
different  local  plasticity  responses  between  the  two  cases.  AEy“ 
and  Ay“3x  values  are  as  large  as  about  70%  in  the  grain  boundary 
region  in  Fig.  6(b)  and  (c).  However,  a  similar  trend  was  also  shown 
for  Actm,  AEy“  and  Ay^^^  between  the  hex-meshed  smooth  and 
the  tet-meshed  smooth  cases  (Fig.  6(d)-(f)),  even  though  Actm. 
A  Ey“  and  A  y“  values  are  much  smaller  for  this  comparison.  This 
indicates  that  the  deformation  constraint  and  associated  deforma¬ 
tion  incompatibility  near  the  grain  boundary  lead  to  different  local 
plasticity  responses  between  the  hex-meshed  smooth  and  the  tet- 
meshed  smooth  grains.  Here,  the  AEy“  and  Ay“3x  values  are  as 
large  as  10%  in  the  grain  boundary  region  (Fig.  6(e)  and  (f)). 

Second,  Actm.  AEy"  and  Ay^^x  values  are  sensitive  to  the 
plastic  properties  of  the  grain.  For  all  comparison  cases,  AEy“ 
and  Ay“3x  values  are  larger  in  grain  2  (multi-slip  oriented  hard 
grain  when  pulled  along  the  tensile  direction)  than  those  in  grain 
3  (single  slip  oriented  soft  grain  when  pulled  along  the  tensile 
direction),  while  Actm  values  are  larger  in  grain  3  than  those  in 
grain  2.  When  compared  to  the  hard  grain,  the  soft  grain  tends  to 
accommodate  local  deformation  gradients  in  an  efficient  fashion 
that  minimizes  AEy“  and  Ay“3x  between  stepped  and  smooth 
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Fig.  6.  Average  percentage  differences  in  local  ctm.  Sy“,  values  between  the  hex-meshed  stepped  and  hex-meshed  smooth  grains  2  and  3  ((a)-(c))  and  between  the 
hex-meshed  smooth  and  tet-meshed  smooth  grains  2  and  3  (d-f)  as  a  function  of  the  normalized  distance  from  the  grain  boundary  at  0.6%  strain. 


grains  and  between  hex-meshed  and  tet-meshed  grains.  In  the  hard 
grain,  however,  the  local  deformation  incompatibility  associated 
with  stair-case  boundaries  seems  to  be  accommodated  by  inducing 
intense  plasticity  activity  in  the  local  grain  boundary  region.  This  is 
probably  why  AEy“  and  values  are  very  large  (particularly, 

larger  E)/“  and  values  in  the  stepped  grain  boundary  region, 
as  indicated  by  solid  circles  in  Fig.  6(b)  and  (c))  in  grain  2.  In  Fig.  6, 
however,  Actm  values  show  the  opposite  trend  -  larger  in  the  soft 
grain  (grain  3)  than  in  the  hard  grain  (grain  2  in  Fig.  6(a)  and  (d)). 
The  elastic  anisotropy  {A  =  2.5  for  grains  2  and  3)  and  associated 
elastic  interactions  was  initially  suspected  to  play  a  role  in  deter¬ 
mining  A(Tm  values  for  the  soft  and  hard  grains  -  for  reference,  the 
calculated  elastic  moduli  for  grains  2  and  3  are  305  and  234  GPa, 
respectively.  However,  in  Fig.  7  Actm  values  are  compared  in  the 
fully  elastic  regime  (at  0.1 5%  strain),  and  these  show  no  significant 
difference  between  grains  2  and  3.  This  means  that  the  Actm  differ¬ 
ence  between  hard  and  soft  grains  in  Fig.  6(a)  and  (d)  arises  in  the 


plastic  regime.  For  the  soft  grain  (grain  3)  in  Fig.  6,  Actm.  AEy“  and 
Ay“ax  values  are  in  the  same  range  (up  to  about  10%  for  the  hex- 
meshed  stepped  and  hex-meshed  smooth  grain  comparison  and 
up  to  about  5%  for  the  hex-meshed  and  tet-meshed  smooth  grain 
comparison),  implying  that  AE/“  and  Ay“ax  tie  directly  to  Actm- 
For  the  hard  grain  (grain  2),  however,  Actm  values  are  far  smaller 
than  AEy“  and  A)/“ax  values  (up  to  about  70%  AE/“  and  A}/"^^ 
and  about  2%  Actm  for  the  hex-meshed  stepped  and  hex-meshed 
smooth  grain  comparison,  and  up  to  about  10%  AE}/“  and  A/"_,x 
and  about  1%  Actm  for  the  hex-meshed  and  tet-meshed  smooth 
grain  comparison).  This  means  that  local  ctm  responses  (perhaps, 
related  local  stress  state  responses)  in  the  hard  grain  seem  to  be  far 
less  sensitive  to  the  grain  boundary  morphology  and  mesh  type, 
compared  to  the  E/“  and  sensitivities  to  the  grain  boundary 
morphology  and  mesh  type  in  the  hard  grain.  This  actually  results  in 
smaller  Actm  values  in  the  hard  grain  than  in  the  soft  grain  (Fig.  6(a) 
and  (d)). 
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Fig.  7.  Average  percentage  differences  in  local  cm  values  (in  the  elastic  regime)  (a)  between  the  hex-meshed  stepped  and  hex-meshed  smooth  grains  2  and  3  and  (b)  between 
the  hex-meshed  smooth  and  tet-meshed  smooth  grains  2  and  3  as  a  function  of  the  normalized  distance  from  the  grain  boundary  at  0.15%  strain. 
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Fig.  8.  Simulated  bicrystal  spheres  visualized  by  (a)  au,  (b)  'Ey"  and  (c)  for  the  hex-meshed  stepped  and  hex-meshed  smooth  and  tet-meshed  smooth  cases. 


As  previously  mentioned,  Actm.  AE)/“  and  A)/“_jx  values  sig¬ 
nificantly  increased  in  the  grain  boundary  region  (as  indicated  by 
broken  and  solid  circles  in  Fig.  6),  compared  to  the  grain  interior. 
All  simulated  bicrystal  spheres  were  visualized  by  their  plasticity 
parameter  values  in  order  to  locate  such  spots  with  large  Actm, 
AE)/“  and  A)/“_jx  values.  Fig.  8  shows  the  grain  surface  and  cross- 
sectional  views  of  ctm  .  and  for  all  simulated  bicrystal  cases 

at  0.6%  strain.  Fig.  8  shows  that  very  large  AE}/“  and  AyJ^^x  values 
in  grain  2  (indicated  by  solid  circles  in  Fig.  6(b)  and  (c))  are  mainly 
distributed  near  the  grain  boundary  region  between  grain  2  and 
grain  1,  as  indicated  by  gray  arrows  in  Fig.  8(b)  and  (c).  Similarly, 
large  Actm  values  in  grain  3  (indicated  by  a  broken  circle  in  Fig.  6(a)) 
were  found  near  the  grain  boundary  region  between  grains  3  and 
1,  as  indicated  by  a  black  arrow  in  Fig.  8(a).  Also,  note  from  Fig.  8 
that  the  overall  global  contours  of  Om.  and  appear  to  be 
similar  between  the  simulation  cases,  even  though  local  gradients 
of  cTivi.  2)/“  and  ^re  somewhat  different  location  by  location 
between  the  simulation  cases. 

4.  Concluding  remarlcs 

In  the  present  study,  parametric  CP-FEM  simulations  were  per¬ 
formed  using  a  bicrystal  (a  soft  grain  and  a  hard  grain)  sphere 
embedded  in  a  cubic  grain  in  order  to  investigate  and  compare 
global  and  local  plasticity  responses  when  the  bicrystal  sphere 
is  represented  by  either  uniform  bexahedral  (cubic)  meshes  with 
stair-case  grain  boundaries,  non-uniform  hexahedral  meshes  with 
smooth  (or  flat)  grain  boundaries  or  tetrahedral  meshes  with 
smooth  (or  flat)  grain  boundaries.  Almost  no  differences  are  found 
in  global  averages  of  ctm.  and  responses  for  all  sim¬ 
ulation  cases  (Fig.  3),  suggesting  that  the  choice  of  the  smooth 
or  stair-case  grain  boundaries  and  the  hexahedral  or  tetrahedral 


meshes  has  a  minimal  influence  on  global  average  responses  at 
the  selected  strains.  However,  the  statistical  distribution  of  ctm. 
E)/“  and  values  showed  that  the  hex-meshed  stepped  case 
causes  more  extreme-value  spots  near  tail  and  head  of  the  dis¬ 
tribution  (Fig.  4).  These  extreme  spots  were  heavily  populated 
in  the  grain  boundary  region  (Fig.  5),  which  is  probably  indica¬ 
tive  of  stair-case  grain  boundaries  as  a  source  and/or  a  sink  for 
excessive  local  extreme  plasticity.  Also,  note  that  the  degree  of 
the  presence  of  those  extreme  spots  in  the  hex-meshed  stepped 
case  depended  upon  the  susceptibility  of  tbe  grain  to  plastic¬ 
ity. 

Larger  AEy“  and  Ay/^^^  values  in  the  hard  grain  (grain  2  in 
Fig.  6(b),  (c),  (e)  and  (f))  indicates  that  the  sensitivities  of  simu¬ 
lated  local  Ey"  and  yl^^x  responses  to  the  choice  of  smooth  or 
stair-case  boundaries  and  hexahedral  or  tetrahedral  meshes  are 
stronger  in  a  hard  grain  than  in  a  soft  grain.  However,  the  soft 
grain  showed  consistent  Ao'm,  AEy“  and  Ay“3x  sensitivities  to 
the  choice  of  the  smooth  or  stair-case  grain  boundaries  and  tbe 
hexahedral  or  tetrahedral  meshes,  while  in  the  hard  grain  the  Activi 
sensitivity  is  significantly  smaller  than  AEy“  and  A/f^^x  sensi¬ 
tivities  to  the  choice  of  the  grain  boundary  representation  and 
the  mesh  type.  Perhaps,  the  further  numerical  study  is  needed 
to  clarify  the  inconsistency  between  the  Actm  sensitivity  and  the 
AEy“  and  Ayl^^^  sensitivities  in  the  hard  grain.  As  expected,  the 
Ao'm,  AE/“  and  Ay"^^  values  between  the  hex-meshed  smooth 
and  the  tet-meshed  smooth  grains  (Fig.  6(d)-(f))  were  fairly  small, 
when  compared  to  those  between  the  hex-meshed  stepped  and 
hex-meshed  smooth  grains  (Fig.  6(a)-(c)).  However,  even  for  the 
former  comparison,  the  simulated  average  local  value  difference 
was  as  large  as  10%  in  Ey“  and  y®^^  and  2%  in  ctivi.  Also,  the  largest 
local  value  difference  was  found  in  the  grain  boundary  region  for 
all  comparisons. 
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Lastly,  note  that  the  current  study  is  limited  to  comparing  plas¬ 
ticity  responses  at  the  relatively  early  stage  of  plasticity  (at  0.6% 
axial  strain),  viz.  the  viscoplastic  regime.  Comparing  Fig.  7  with 
Fig.  6(a)  and  (d)  indicates  that  the  trend  in  Actm  between  differ¬ 
ent  grain  boundary  representations  and  mesh  types  is  different  in 
the  fully  elastic  regime  and  in  the  viscoplastic  regime.  It  would 
be  worthwhile  to  investigate  how  the  trend  in  Ao'm,  AEy“  and 
^/max  between  different  grain  boundary  representations  and  the 
mesh  types  changes  at  a  large  plastic  strain,  which  is  the  case  for 
metal  forming  simulations.  Such  a  numerical  study  would  provide 
useful  information  for  evaluating  the  sensitivity  of  local  plasticity 
responses  to  different  choices  of  the  grain  boundary  representa¬ 
tions  and  mesh  types  in  metal  forming  simulations. 
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